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ABSTRACT 


The  objective  of  this  study  is  to  verify  the  Symmetrical 
Number  System  (SNS)  undersampling  receiver  architecture  using 
software  and  investigate  implementation  issues  using  Digital 
Signal  Processing  (DSP)  hardware.  In  the  software  design,  a 
MATLAB  program  is  written  to  determine  a  single  sinusoidal 
input  frequency  using  this  receiver  architecture.  Each 
channel  of  the  SNS  under satr^ling  receiver  consists  of  a  low 
speed  ADC,  a  discrete  Fourier  transform  followed  by  a 
constant  threshold  device  to  detect  the  signal's  frequency 
bin.  The  detected  frequency  bins  are  then  recombined  in  a 
SNS-to-decimal  algorithm  to  recover  the  frequency  of  the 
signal .  Error  rate  performance  in  a  Gaussian  noise 
environment  at  the  input  stage  is  evaluated.  In  the  hardware 
design,  a  sinusoidal  waveform  is  digitized,  discrete  Fourier 
transformed  and  converted  from  the  SNS  format  to  a  decimal 
value  using  a  single  channel  digital  signal  processor. 
Implementation  difficulties  and  design  issues  are  discussed. 
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I .  INTRODUCTION 


A.  UNDERSAHPLIN6 

The  digitization  of  a  signal  is  usually  governed  by  the 
Nyquist  theorem  where  the  sampling  frequency  is  at  least 
twice  the  signal  bandwidth.  The  Nyquist  theorem  however, 
places  a  limitation  only  on  the  information  that  can  be 
derived  from  a  single  set  of  digitized  data  [Ref.  1].  If  the 
sampling  frequency  is  less  than  twice  the  bandwidth  of  the 
signal  being  digitized,  aliasing  and  consequently  ambiguities 
occur.  With  additional  information  however,  ambiguous 
frequency  components  due  to  undersampling  may  be  resolved. 
Such  information  may  come  from,  for  example,  trial  sampling 
periods.  Rader  [Ref.  2]  described  how  trial  sampling  periods 
can  be  used  to  recover  periodic  signals.  The  trial  sampling 
period  which  yields  the  waveform  of  smallest  variation  is 
considered  to  be  the  correct  period  and  the  resulting 
waveform  the  correct  waveform. 

Pace,  Leino  and  Styer  [Ref.  3]  examined  the  relationship 
between  the  Discrete  Fourier  Transform  (DFT)  and  the 
Symmetrical  Number  System  (SNS)  as  a  means  of  resolving 
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single  frequency  undersampling  aliases.  They  showed  that  the 
DFT  encodes  the  frequency  information  of  a  signal  in  a  format 
that  is  in  the  same  form  as  the  SNS.  In  addition,  they 
proved  analytically  that  aliases  resulting  from  undersampling 
a  single- frequency  signal  could  be  resolved  using  2  or  more 
channels.  Each  channel  in  a  SNS  undersartpling  receiver 
contains  a  low  speed  ADC,  a  DFT  and  a  threshold  device  to 
detect  the  input  signal  bin  number  in  the  frequency  domain. 
The  bin  numbers  from  each  channel  are  then  recombined  to 
resolve  the  signal's  frequency. 

B .  PRINCIPAL  CONTRIBUTIONS 

First,  this  thesis  verifies  the  SNS  undersampling  theory 
advanced  by  Pace,  Leino  and  Styer  [Ref.  3] .  An  algorithm  is 
written  and  coded  in  MATLAB  to  prove  the  methodology  and  to 
show  that  the  frequency  of  an  under  sampled  signal  can  be 
accurately  measured.  The  algorithm  is  also  simulated  in  a 
Gaussian  noise  environment.  Error  rates  for  the  different 
noise  levels  are  obtained  as  a  function  of  the  signal  to 
noise  ratio.  Since  the  Fast  Fourier  Transform  (FFT)  is  not 
suitable  for  computing  DFTs  in  this  application,  alternative 
methods  are  suggested  for  real-time  applications. 
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Second,  possible  hardware  implementation  problems  are 
investigated  based  on  a  Digital  Signal  Processing  (DSP) 
platform.  Several  problems  were  encountered:  the  need  for 
stable  sampling  frequencies,  large  memories  and  alternative 
methods  for  computing  DFT  for  fast  response  time. 
Integration  into  future  EW  receivers  must  take  these  factors 
into  consideration. 

Undersampling  offers  several  advantages  [Ref.  4]  .  It 
allows  the  resolution  of  very  high  frequencies  in  EW 
receivers  using  low  speed  ADCs.  This  is  especially  so  if 
several  SNS  channels  are  used.  In  particular,  the  use  of 
undersampling  in  the  design  of  receivers  will  reduce  their 
cost  and  complexity. 

C.  THESIS  ORGANIZATION 

In  Chapter  II,  the  relationship  between  the  SNS  and  the 
digital  frequency  domain  as  mapped  by  the  DFT  is  examined  as 
a  means  of  resolving  single- frequency  undersattqpling 
ambiguities.  It  shows  how  the  frequency  of  a  signal  that  is 
undersampled  at  two  different  sampling  frequencies  (two- 
channel)  can  be  determined.  In  order  to  use  lower  sampling 
frequencies,  the  two- channel  case  can  be  extended  to  three  or 
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more  channels.  In  particular  the  three-channel  case  is 
discussed. 

In  Chapter  III,  algorithms  for  the  two-channel  and 
three-channel  receivers  are  developed  and  coded  in  MATLAB  to 
measure  the  frequency  of  an  incoming  signal.  Each  section  of 
the  software  is  explained  in  detail.  Results  are  obtained 
based  on  different  Gaussian  noise  levels. 

A  feasibility  study/design  for  the  two-channel  case  is 
carried  out  in  Chapter  IV  using  a  DSP  development  kit.  The 
suitability  of  using  a  DSP  platform  and  its  associated 
problems  are  discussed. 

Chapter  V  states  some  conclusions  and  recommendations 
for  future  research. 
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II .  BACKGROXniD  INFORMATION 


A.  INTRODUCTION 

Digitization  of  a  signal  is  usually  governed  by  the 
Nyquist  criterion  when  the  input  signal  is  bandlimited  to  0  < 
f  <  fs/2  where  is  the  sampling  frequency.  For  higher 
frequencies  (i.e.  f  >  t^/2) ,  the  process  of  undersampling 
gives  rise  to  ambiguities.  However,  with  additional 
information  (or  channels) ,  the  frequency  components  f  >  tj2 
can  be  resolved. 

Pace,  Ramamoorthy  and  Styer  [Ref.  5]  showed  that  the 
discrete  Fourier  transform  (DFT)  naturally  encodes  the 
frequency  information  of  a  signal  in  the  same  format  as  the 
symmetrical  number  system  (SNS) .  Consequently,  aliases  from 
undersampling  can  be  resolved  using  this  method.  The  theory 
set  forth  is  elaborated  in  [Ref.  3] . 

B.  DISCRETE  FOURIER  TRANSFORM  (DFT) 

Since  all  signals  consist  of  sinusoids,  for  simplicity, 
a  single  frequency  sinusoidal  waveform  is  used  for  analysis. 
Assume  the  sinusoidal  signal  is 

x(t)  =  2  cos  ®  t 
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(1) 


and  after  satrpling 

x(n)  =  2  cos  ®  n. 


The  DFT  of  x{n)  is  given  by  [Ref.  6] : 


(2) 


N-l 


X(k)  =  k  =  0,1,. ..,N  - 1. 


n=0 


(3) 


Applying  the  DFT  to  x(n)  results  in  a  discrete  spectrum  where 
|X(k)  1^  is  the  energy  contained  in  the  signal  at  each  digital 
frequency  ®=2xk/N.  The  spectrum  X(k)  has  N  indices  with  the 
digital  frequency  of  each  index  given  by: 


for  N  even 


and 


(4) 


0,271— ,...,271  — 

N  N 


N 


N 


N 


for  N  odd. 


(5) 
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The  analog  frequency  corresponding  to  each  index  is  obtained 
by  multiplying  each  value  by  f^.  Since  signals  with  digital 
frequencies  in  the  range  n  <  &  <  2%  are  indistinguishable 
from  signals  with  digital  frequencies  0  <  ©  <  ti,  the  digital 
frequency  of  each  index  can  also  be  written  as: 


n,,  '  (N/2-1)  -1,1 

0,271  ,...,271  ,27c  ,,..,271  ,27C 

N  N  N  N  N 


for  N  even 


and 


(6) 


1  ^  |N/2j  ^  Ln/2J  ^  2  ^  r 
N  N  N  N  N 


for  N  odd. 


(7) 


where  LxJ  is  the  floor  function  and  represents  the  greatest 
integer  less  than  or  equal  to  x.  Thus  the  spectrum  X(k) 
resolves  into  N  integer  indices  and  incoming  signals  will  map 
into  unique  bins: 


N  N_ 

2  ’  2  ’ 


for  N  even. 


(8) 


N 

N 

0,1,..., 

7; 

_2_ 

,...,2,1 

for  N  odd. 


(9) 
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For  example,  for  N  =  5  (£3  =  5  Hz  and  the  sampling  duration 
Ti  is  1  second)  ,  the  output  bins  after  the  DFT  are  [0122 
1]  for  input  frequencies  of  [0  1  2  3  4]  Hz.  These  DFT  bins 
are  repeated  for  higher  frequencies  as  illustrated  in  Figure 
1.  In  this  figure  the  abscissa  corresponds  to  the  incoming 
frequency  and  the  ordinate  corresponds  to  the  bin  into  which 
the  signal  is  resolved. 
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Figure  1 :  DFT  bin  mapping  for  input  frequencies  f = 

0  to  10  for  N  =  5  (f^  =  5  Hz  sampling  for  1  second) . 
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C.  THE  SY21METRICAL  NUMBER  SYSTEM  (SNS) 

The  SNS  is  coitposed  of  a  number  of  pairwise  relatively 
prime  (PRP)  moduli.  The  integers  within  each  SNS  modulus 
however,  are  derived  from  a  symmetrically  folded  waveform. 
The  symmetrically  folded  waveform  corresponding  to  each  SNS 
PRP  moduli  (mj  ,  has  a  folding  period  equal  to  the  modulus . 
The  integer  values  within  each  SNS  modulus  are  derived  from  a 
tiid- level  quantization  of  the  symmetrical  folding  waveform. 
The  formal  definition  of  a  symmetrical  residue  is  given 
below: 

Definition:  For  an  integer  h  such  that  0  <  h  <  m 


=min  (h,  m-h} 

(10) 

If  this  function  is  extended  periodically  with  period  m, 
that  is, 

^h  +  nm  “ 

(11) 

where  n  e  {0,±1,±2, ...}  then  x^  is  called  a  symmetrical 
residue  of  (h+nm)  modulo  m.  For  m  even,  let  x  be  the  row 
vector 
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For  m  odd,  let  x  be  the  row  vector 


"" 

m 

m 

x  = 

0,1,..., 

7; 

_Y_ 

,...,2,1 

(13) 


where  LxJ  again  represents  the  floor  function  resulting 
in  the  greatest  integer  less  than  or  equal  to  x.  These 
two  vectors  consist  of  the  symmetrical  remainder 
elements  x^,  0  <  h  <  m. 


D.  RELATIONSHIP  BETWEEN  DFT  AND  SNS 

From  the  aibove,  it  is  obvious  that  the  DFT  maps  real 

signals  naturally  into  the  SNS.  That  is,  in  Section  C,  if  we 
let  the  modulus  m  represent  the  san^ling  frequency  multiplied 
by  the  saitpling  time  (i.e.,  f^Ti)  ,  then  equations  (12)  and 
(13)  are  in  the  same  form  as  equations  (8)  and  (9)  where  N= 
fgTi.  Thus  the  SNS  provides  a  convenient  framework  for 
undersatrping  signal  analysis. 


Table  1  displays  the  input  frequencies  and  the  resulting 
DFT  bins  for  sampling  frequencies  5  Hz  and  6  Hz  respectively. 


Input  Frequency 

DFT 

Bins 

f 

f^=  5  Hz 

f^=  is  Hz  I 

0 

} _  ^ 

0 

0  ! 

1 

1 

1  1 

2 

2 

2  ! 

r  3 

2 

3  i 

r~  4 

1  1 

2  i 

1  5 

i -  .  _ 

1  0  j 

1  i 

i  6  1 

1  1 

0  1 

1  7  i 

:  2  1 

1  i 

1  8  I 

■  ■  r*  ■ 

2  j 

2  1 

Table  1:  Input  Frequency  and  Resulting  DFT  Bins  for  2 
Channel  Example . 


The  frequencies  are  resolved  as  described  in  equations  (12) 
and  (13)  .  By  considering  two  or  more  channels,  it  is 
possible  to  unambiguously  resolve  the  signal  frequencies  in 
the  dynamic  range  determined  by  the  SNS.  One  method  is  to 
devise  a  look-up  table  similar  to  that  shown  in  Table  1 . 
However  this  method  is  inefficient  for  high  frequencies; 
large  memories  are  required.  An  alternative  method  is 
described  below: 

Suppose  there  are  r  channels  and  the  incoming  frequency 
is  within  the  dynamic  range  of  the  system.  To  carry  out  the 
SNS-to-decimal  conversion,  we  need  to  solve  f  =  ai(mod  iiii) 
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for  i 


where  Hi  is  the  corresponding  detected  DFT 


bin  for  each  nii.  The  Chinese  Remainder  Theorem  states  that 
there  is  a  unique  solution  modulo  M=m3L*m2  ...*mj,.  A  standard 
method  of  solution  is  to  find  integers  such  that  M*bi/mi  = 
l(mod  trii)  where  i  =  l,2,....r  in  which  case  the  solution  is  f 
=  M*bi*ai/mi  +  M*b2*a2/m2  +  .  . .  +  iyi*br*aj./irv  (mod  M)  .  In  Sections 
F  and  G  below,  examples  are  given  to  illustrate  this 
calculation. 

E.  D1!NAMIC  RANGE  OF  THE  SHS 

Let  ttij,  .  .  . ,  m^  be  r  pairwise  relatively  prime  moduli, 

then  the  dynamic  range,  D  (0:D-1)  of  a  SNS  system  is  given  as 
follows : 

•  If  all  the  moduli  are  odd,  then  the  dynamic  range  of 
the  system  is 


D  =  min 

[  ^  1=1  ^  I=j+1 

(14) 

where  j  ranges  from  1  to  r-1  and  range  over 

all  permutations  of  {l,  2 , 3 , r} .  For  example,  for  a 
two-channel  case  with  m^  =  5,  ttij  =  7, 
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^  .  fm,  m, 

D  =  mim  — L-) — 2 

1  2  2  . 


or  D  =  6. 


For  a  three-channel  case  with  tni=  3,  in2  =  5,  103  =  7, 


D--mm{m,  +m2m3,m2  +m,m3,mj  +mim2} 


or  D=22. 

of  the  moduli  (mj^)  is  even,  then  the  dynamic 
range  of  the  system  is 


D  =  min|l^I^m,,+n 


m. 


2  i=j+i 


(15) 

where  j  ranges  from  1  to  r-l  and  ...m;  range  over 
all  permutations  of  {2,3,...,r}.  For  example,  for  a 
two- channel  case  with  m^  =  6  nij  =  5, 


.  f 


m, 


D  =  min|-^+m 


or  D  =  8, 
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For  a  three-channel  case  with  tni=  8,  tn2  =  5,  103  =  7, 


m,  m,  nij 

D  =  min<^ - — m,  +  m,,  — m,  +  m 

[2  "  2  "  2  ^ 

or  D=27. 

Clearly,  the  dynamic  range  of  an  SNS  system  with  one 
even  modulus  is  superior  to  that  using  all  odd  moduli . 
Moreover,  the  greater  the  number  of  channels,  the  greater  the 
dynamic  range. 

F.  THE  TWO-CHAHNEL  CASE 

Figure  2  shows  the  block  diagram  of  a  two -channel 
receiver  architecture  to  determine  a  single  frequency  f.  In 
this  architecture  the  ADC  sampling  frequencies  f^^  and  fgj  are 
relatively  prime  and  =  1.  The  DFT  outputs  are  thresholded 
to  detect  the  frequency  bins  of  the  signal.  The  detected 
frequency  bins  a^  and  aj  are  then  used  by  the  SNS -to -decimal 
algorithm  to  determine  the  frequency  of  the  input  signal. 
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frequency  bins  and  a^  are  then  used  by  the  SNS-to-decimal 
algorithm  to  determine  the  frequency  of  the  input  signal. 


f 


Figure  2:  Block  Diagram  of  a  Two  Channel  Receiver 
Architecture . 


Let  till  =  f^i  and  ttij  =  f^j  and  suppose  that  the  incoming 
frequency  is  within  the  dynamic  range  of  the  system.  From 
Section  D,  we  need  to  solve  f  s  ai(mod  m^)  and  f  s  a^  (mod  la,)  . 
The  two  congruence  equations,  f  s  ai(mod  m^)  and  f  =  a^  (mod 
mi)  are  solvable  only  if  the  greatest  common  divisor  of  irii 
and  mi  divides  (ai  -  a^)  ,  a  generalization  of  the  Chinese 
Remainder  Theorem  [Ref.  7].  To  solve  for  f,  the  diophantine 
equation 


p*mi+q*mi  =  (ai  -  aj 


(16) 
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must  be  solved  for  p  and  f  is  then  calculated  from  the 
ec[uation 

f  =  ai+p*mi.  (17) 

The  code  for  this  algorithm  is  shovm  in  Appendix  A. 

For  example,  for  sartqpling  frequencies  5  and  6,  m^  and  it^ 
have  values  of  5  and  6  respectively  (Tj  =  1)  .  If  the  signal 
is  resolved  into  bins  a^  (=2)  and  aj  (=1)  after  the  DFT,  p 
is  found  to  have  a  value  of  1  and  q  is  found  have  a  value  of 
-1.  Thus,  the  input  frequency  from  (17)  is  2+1*5  =  7.  This 
can  also  be  verified  as  shown  in  Table  1. 

G.  THE  THREE-CHAHNEL  CASE 

Figure  3  shows  the  block  diagram  of  a  three- channel 
receiver  architecture  to  determine  a  single  frequency  f. 
Similar  to  the  two- channel  case,  the  ADC  sampling  frequencies 
fai»  ^82/  fs3  pairwise  relatively  prime  and  T^  =  1.  The 

DFT  outputs  are  thresholded  to  detect  the  frequency  bins  of 
the  signal.  The  frequency  bins  a^,  and  are  then  used  by 
the  SNS-to-decimal  algorithm  to  determine  the  frequency  of 
the  input  signal. 
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Figure  3:  Block  Diagram  of  a  Three  Channel 
Receiver  Architecture . 


In  the  three-channel  solution,  let  and 

=  fs3  and  suppose  that  the  incoming  frequency  is  within  the 
dynamic  range  of  the  system.  We  need  to  solve  f  s  a^Cmod  nii) 
and  f  =  aj  (mod  tiij)  and  f  =  aj  (mod  trig)  .  Using  the  Chinese 
Remainder  Theorem  and  the  Euclidean  algorithm,  the  method  of 
solution  is  to  find  integers  bi  such  that  M*bi/mi  =  l(mod  111^) 
where  i  =  1,2,  and  3  and  M  =  mj^  *  m^  *  m3 .  The  solution  is 
then  f  =  ±  M*bi*aiM  ±  M*b3*a3/m3  ±  M*b3*a3/m3  (mod  M)  where  f 
is  the  frequency  which  falls  within  the  dynamic  range  D  of 
the  system. 
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For  example,  let  =  5,  nij  =  6  and  ti^  =  7,  so  that  M  = 
210  and  D  =  22 .  Suppose  that  the  signal  is  resolved  into 
bins  a^  (=  1)  ,  a.2  (=  2)  and  (=2)  after  the  DFT.  For  the 
three- channel  case  the  bi  values  must  be  found.  Here,  b^^  b2 
and  bs  are  found  to  be  -2,  -1,  and  -3  respectively.  Thus  f  = 
±  210(-2)(l)/5  ±  210(-l)(2)/6  +  210(-3)(2)/7  mod(210)  and  we 
must  choose  the  solution  that  falls  within  the  SNS  dynamic 
range  D  =  22  [0:21].  The  correct  combination  f  s  84  -70  + 

180  mod  (210)  =  194  mod  (210)  .  Although  194  is  out  of  the 

dynamic  range,  210  -  194  =  16  is  in  the  dynamic  range  so  that 
f  =  16  is  the  correct  frequency. 


H.  NOISE  CONSIDERATIONS 

For  a  sinusoidal  waveform,  the  Signal  to  Noise  Ratio 
(SNR)  is  defined  as 


p 

SNR  =  -V 
2c^ 


(18) 


where  P  is  the  power  of  the  signal  and  is  the  noise  power. 
Assuming  a  signal  power  of  one,  the  noise  power  and  amplitude 
are  given  by 
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(20) 


This  a  is  multiplied  by  a  normally  distributed  random  number 
sequence  of  zero  mean  and  unit  variance  and  added  to  the 
input  signal  as  noise.  The  simulation  results  are  given  in 
Chapter  III. 
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III.  SOFTWARE  DESIGN  AND  RESULTS 


A.  TWO-CHANNEL  ALGORITHM 

The  two- channel  case  was  described  in  Chapter  II.  An 
algorithm  was  constructed  based  on  Figure  2 .  The  software 
given  in  Appendix  A  can  be  divided  into  the  following 
sections : 

•  Initialization.  This  section  obtains  all  the 
parameters  (number  of  iterations,  input  frequency, 
sampling  frequencies,  quantization  levels)  required. 

•  Iteration  loop.  This  section  consists  of  a  loop 
(with  an  initial  count  of  zero)  to  count  the  number 
of  errors . 

•  Creation  of  Waveform.  Based  on  the  input  frequency, 
a  sinusoidal  waveform  is  created  with  noise  added. 

•  Sampling  and  Quantization.  The  waveform  is  then 
saitpled  at  two  different  frequencies  and  quantized 
using  a  14 -bit  ADC. 

•  Windowing.  A  rectangular  window  operation  of  width  N 
=  fg  *  Ti  =  fg  (the  total  sampling/integration  time  is 
taken  to  be  one) is  carried  out. 
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•  DFT  Operation.  A  DFT  is  then  carried  out  on  each 
sample,  taking  only  the  first  half  of  the  DFT  output. 
The  formula  used  for  the  DFT  process  is  a  simple  pair 
of  nested  loops. 

•  Bin  Detection.  A  non-adaptive  (constant)  threshold 
bin  detector  is  then  used  to  find  the  bin  with  the 
maximum  value  for  each  DFT  output. 

•  SNS-to-Decimal  Algorithm.  The  SNS-to-decimal 
algorithm  as  described  in  Chapter  ll  is  then  used  to 
calculate  the  incoming  frequency. 

A  flow  diagram  of  this  algorithm  is  illustrated  in 

Figure  4 .  The  MATLAB  code  can  be  found  in  Appendix  A . 
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Figure  4 :  Two  Channel  Algorithm 
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B.  TESTING  OF  TWO-CHANNEL  SYSTEM 


To  test  the  two- channel  case  (sinusoidal  signal  without 
noise)  ,  the  program  is  run  with  the  following  input  and 
sampling  frequencies  shown  in  Table  2. 


f 

fsl 

fs. 

Dynamic  Range 

Remarks 

7 

5 

8 

0:8 

Low  input  frequency 

100 

97 

98 

0:145 

Consecutive  sampling 
frequencies 

1040 

547 

1200 

0:1146 

Sampling  frequencies 
far  apart 

12125 

12671 

12919 

0:12794 

High  input  frecjuency 

Table  2 :  Tested  Input  and  Sanpling  Frequencies 


For  exanple,  with  input  signal  frequency  at  7  Hz  as 
shown  in  Figure  5,  the  sampled  signals  at  5  Hz  and  at  8  Hz 
are  shown  in  Figures  6  and  7  respectively.  The  DPT  output 
for  the  two  saitples  are  shown  in  Figures  8  and  9 .  The 
resultant  bins  of  the  first  halves  of  Figure  8  and  9  are  then 
supplied  to  the  SNS-to-decimal  algorithm  to  be  converted  to 
the  input  frequency  of  7  Hz. 
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Figure  8:  DFT  output  with  fsl=5  Hz. 
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Figure  9:  DFT  output  with  fs2=8  Hz 


It  is  found  that  if  one  of  the  sampling  frequencies  was 
the  same  as  the  input  frequency,  the  algorithm  failed.  This 
is  because  the  resulting  samples  due  to  the  same  sampling 
frequency  will  consist  of  zeros.  This  problem  can  be  solved 
by  using  at  least  two  sets  of  sampling  frequencies .  Apart 
from  this,  the  algorithm  works  well  in  this  noise-free  (high 
signal -to-noise  ratio)  environment. 

C.  SIZSDLATION  PARAMETERS  FOR  TWO-CHAMMEL  CASE 

To  obtain  the  error  rates  in  a  noisy  environment,  the 
two-channel  software  is  run  with  the  following  parameters: 

•  Number  of  iterations,  num  =  10000 

•  Signal  to  Noise  Ratio,  SNRDB  =  -30  to  30  dB 

•  ADC  resolution,  bit  =  14 

•  Input  and  sampling  frequencies  as  shown  in  Table  3. 


f 

fsx 

9 

10 

11 

90 

91 

92 

900 

901 

902 

9000 

9001 

9002 

Table  3:  Input  and  Sampling  Frequencies. 

D.  RESULTS  FOR  TWO- CHANNEL  CASE 

The  results  obtained  are  shown  in  Figure  10 . 
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Figure  10 :  Error  Rates  vs .  SNR  for  two-channel  system 
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The  following  observations  are  made: 


•  As  expected,  the  error  rates  improve  as  the  SNR 
increases.  A  tradeoff  between  SNR  and  error  rate  is 
required. 

•  Improvements  in  error  rates  were  obtained  when  higher 
frequencies  were  used.  This  is  because  at  higher 
frequencies,  higher  satipling  frequencies  are 
rec[uired.  This  leads  to  a  higher  N-point  DFT  (higher 
gain)  which  is  less  affected  by  noise . 

•  However,  at  higher  frequencies,  the  time  taken  to 
compute  the  DFT  was  much  longer.  To  reduce  the  time 
taken,  the  following  methods  can  be  in^lemented: 

•  If  N  is  highly  composite  (factorable  into  powers  of 
many  small  prime  factors,  preferably  primes  <  10) , 
use  a  "mixed-radix"  FFT  implementation. 

•  If  N  is  prime,  or  contains  very  large  prime 
factors,  use  the  "chirp-z"  transform. 

•  Use  three  or  more  channels  in  the  receiver.  A 
three-channel  receiver  has  a  higher  dynamic  range 
for  the  same  magnitude  of  sampling  frequencies . 
For  example,  a  two-channel  receiver  with  sampling 
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frequencies  6  and  7  has  a  dynamic  range  of  [0:9] 
while  a  three-channel  receiver  with  sampling 
frequencies  of  5,  6  and  7  has  a  dynamic  range  of 

[0:21]  . 

E.  THREE -CHAEMEL  ALGORITHM 

The  three- channel  algorithm  is  similar  to  the  two- 
channel  algorithm  as  shown  in  Figure  11.  The  MATLAB  code 
can  be  foxmd  in  Appendix  A. 
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Figure  11:  Three  Channel  Algorithm 
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F.  TESTING  OF  THREE  CHANNEL  ALGORITHM 

To  test  the  three- channel  case,  the  program  is  run  with 
some  of  the  following  input  and  sampling  frequencies  in  Table 

4. 


f 

fsl 

fs. 

f.3 

Dynamic 

Range 

Remarks 

13 

5 

6 

7 

0:21 

Low  input  frequency 

100 

17 

18 

19 

0:171 

Consecutive 

sampling 

frequencies 

1040 

17 

91 

919 

0:1232 

Sampling 
frequencies  far 
apart 

12125 

90 

929 

937 

0:42741 

High  input 
frequency 

Table  4:  Tested  Input  and  Sampling  Frequencies. 


Apart  from  the  anomaly  discussed  in  the  two- channel 
case,  the  algorithm  works  well  in  this  noise-free  (high 
signal-to-noise  ratio)  environment. 

G.  SIMULATION  PARAMETERS  FOR  THREE-CHANNEL  CASE 

To  obtain  the  error  rates  in  a  noisy  environment,  the 
three-channel  software  is  run  with  the  following  parameters: 

•  Number  of  iterations,  num  =  10000 
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Signal  to  Noise  Ratio,  SNRDB  =  -30  to  30  dB 


•  ADC  resolution,  bit  =  14 

•  Input  and  sattpling  frequencies  as  shown  in  Table  5. 


f 

fel 

fs. 

fs3 

9 

5 

7 

11 

90 

13 

14 

17 

900 

41 

42 

43 

9000 

141 

142 

143 

Table  5:  Input  and  Sampling  frequencies. 

H.  RESULTS  FOR  THREE-CHANNEL  CASE 

The  results  obtained  are  shown  in  Figure  12 . 
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Figure  12 :  Error  Rates 


for  three-channel  system 


Comparing  the  two-channel  and  three-channel  cases,  the 
following  observations  can  be  made: 

•  The  three- channel  system  is  much  faster  than  the  two- 
channel  system  since  the  DFTs  required  are  smaller 
due  to  the  smaller  sampling  frequencies. 

•  However  the  results  for  the  two- channel  system  with 
noise  are  better.  For  example  to  achieve  a 
relatively  error-free  system  for  a  frequency  of  9000 
Hz,  the  two-channel  case  requires  only  -22  dB. 
However,  the  three-channel  case  requires  at  least  -4 
dB. 
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IV.  HARDWARE  DESIGN  AND  FINDINGS 


A.  INTRODUCTION 

In  the  last  chapter,  the  instantaneous  measurement  of 
frequency  using  the  SNS-to-decimal  algorithm  was  verified. 
There  is  a  need  to  investigate  the  implementation  of  the 
algorithm  in  hardware.  Digital  Signal  Processing  (DSP) 
hardware  was  selected  for  the  following  reasons: 

•  A  major  part  of  the  algorithm  is  the  processing  of 
DFTs  which  is  a  digital  signal  processing  task  well 
suited  to  be  carried  out  by  DSP  hardware. 

•  DSP  hardware  provides  a  fast  way  to  inclement  the 
algorithm.  The  DSP  development  kit  is  easy  to  learn, 
program  and  simulate.  It  is  ideal  for  this 
application  to  investigate  hardware  problems  and 
limitations . 

•  Cost  consideration:  the  development  kit  plus  tools 
cost  $1500; 


•  EW  receivers  are  likely  to  incorporate  DSP  hardware. 


B.  TI  TMS320C54X  DSP  DEVELOPMENT  KIT 

The  TMS32054C54X  DSKplus  [Ref.  8-12]  is  a  low  cost  DSP 
Starter  kit  that  gives  a  designer  a  working  knowledge  of  DSP 
code  to  build  DSP  based  systems .  The  development  kit 
contains  a  stand-alone  application  board  that  can  be 
connected  to  the  PC.  It  executes  code  in  real  time  at  40 
MIPS  while  the  Windows-based  debugger  analyzes  it  line-by- 
line,  displaying  internal  DSP  register  information  in 
multiple  windows  and  in  real  time.  It  has  an  Analog 
Interface  Circuit  for  the  input  of  signals.  The  board's 
communication  interface  enables  the  creation  of  C54x  DSP  code 
and  host  PC  code.  Moreover,  the  hardware  enables  the  use  of 
expansion  slots  for  adding  memory,  peripherals  such  as 
interface  logic,  other  DSPs  etc.  The  developed  code  can 
eventually  be  loaded  into  a  resident  DSP  processor,  which  may 
be  part  of  a  EW  receiver  architecture.  Figure  13  shows  a 
block  diagram  of  the  development  kit.  A  more  detailed 
description  of  the  kit  can  be  found  in  Appendix  B. 
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Parallel  Port  to  PC  To  Other  Interface  Ports 

Figure  13:  Block  Diagram  of  DSP  Hardware. 


C.  SOFTWARE 

The  software  for  the  two- channel  case  described  in 
Chapter  II  (Figure  2)  is  written  using  the  DSP  development 
kit.  The  software  (found  in  Appendix  C)  is  coded  in  'C' 
language/assembly  language  and  converted  to  the  C54x  assembly 
language  (if  required)  prior  to  execution: 

•  Firstappl . c/Firstapp2 . c .  These  two  programs  poll  the 
input  channel  and  sait^le  the  input  signal  at  the  two 
sampling  frequencies  respectively. 

•  Hostappl.cpp/Hostapp2 .cpp.  These  two  programs 

display  the  samples  of  the  signals  based  on  the  two 
sampling  frequencies  and  save  the  data  in  text  files. 
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•  Main.c.  This  program  reads  the  data,  executes  the 
DFT,  obtains  the  largest  values  for  the  two  channels 
and  then  carries  out  a  SNS-to-decimal  conversion. 
These  programs  were  run  individually  and  consecutively. 

D.  TESTING  AND  RESULTS 

Using  data  generated  by  MATLAB,  the  main  program  was 
tested  successfully  in  the  development  kit.  The  programs 
were  then  run  with  an  input  frequency  of  126  Hz  and  satrpling 
frequencies,  125  Hz  and  128  Hz.  Results  obtained  were 
intermittent  i.e.,  correct  results  were  not  always  obtained. 
A  frequency  counter  and  an  oscilloscope  were  set  up  and  it 
was  found  that  the  sampling  frequencies  were  not  stable. 
Testing  with  different  frequencies  did  not  inprove  the 
results . 

E .  PROBLEMS 

Several  problems  were  encountered  during  the 

investigation : 

•  Stability  of  Satr^jling  Frequencies.  The  development 
kit  carries  out  frequency  division  of  the  master 
oscillator  to  obtain  the  sampling  frequencies. 
Unfortunately,  the  crystal  oscillator  has  a 
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resolution  of  5-10  Hz.  This  is  unacceptable  as  a 
shift  of  1  Hz  in  the  sampling  frequency  will  cause 
erroneous  results.  Moreover,  the  fact  that  the 
sampling  frequencies  are  factors  of  the  oscillator 
frequency  and  that  they  need  to  be  pairwise 
relatively  prime  severely  limits  the  choice  of 
frequencies.  A  possible  solution  is  to  obtain  the 
sampling  frequencies  directly  from  stable  signal 
sources . 

•  DFT.  For  higher  frequencies,  the  execution  of  the 
DFT  takes  a  long  time.  Several  solutions  were 
suggested  and  discussed  in  the  previous  chapter. 

•  Memories.  Insufficient  memory  error  messages  were 
encountered  when  high  frequencies  were  used.  The 
same  messages  occurred  when  attempts  were  made  to  run 
the  routines  together.  More  memories  and/or  more 
efficient  DFT  algorithms  are  required. 
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V.  CONCLUDING  REMARKS 


The  main  contribution  of  this  thesis  is  the  verification 
of  the  relationship  of  the  DFT  to  the  SNS  to  resolve 
undersan^ling  ambiguities  and  the  investigation  of  hardware 
irr5>lementation  issues  using  a  DSP  platform.  Error  rates  for 
different  SNR  are  also  obtained. 

The  use  of  undersampling  technique  using  the  SNS  to 
measure  frequency  is  a  viable  method  to  inclement  in  a  EW 
receiver  architecture.  However,  the  need  for  faster  DFT 
computation  and  stable  saitpling  frequencies  must  be  taken 
into  account  before  they  can  be  considered  for  incorporation 
into  EW  receivers.  There  is  also  a  trade-off  between  the 
number  of  channels  and  SNR.  For  faster  response,  a  multi¬ 
channel  case  is  recommended;  but  a  higher  SNR  is  required. 
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APPENDIX  A 


MATLAB  CODE  FOR  SOFTWARE  ALGORITHM 


%  Thesis  Project 
% 

%  Two  Channel  Receiver 
% 

%  Note:  The  sampling  frequencies  should  be  relatively  prime 


clear  all; 


%  Initialization 


num=input ( *  Enter  Number  of  iterations : ' ) ; 
f =input { *  Enter  Input  Frequency: ' ) ; 
fsl=input (* Enter  Sampling  Frequency  1:*); 
fs2=input (* Enter  Sampling  Frequency  2:*); 
fpl=fopen( ’c:\matlab\bin\thesis\result.dat  * 


%  Number  of  iterations 
%  Frequency  of  signal 
%  Sampling  frequency  1 
%  Sampling  frequency  2 
*at');  %  Store  results 


%  Quantization  levels 


bit=14; 

qnlevel=2  "^bi  t  - 1  ; 
q=2/qnlevel; 


%  No.  of  quantization  levels 
%  quantization  size 


for  SNRDB=-30:2:30 


%  Set  Signal  to  Noise  Ratio 
%  from  -30  dB  to  30  dB 


count =0; 


%  Error  Count 


for 


i=l :num 

SNR=10^(SNRDB/10) ; 
sigmasq=l/2/SNR; 


%  Convert  to  non-dB  units 
%  Noise  normalization  assuming 
%  signal  power  of  1 


t=(0: .001:1) ; 
sig=sin(2*pi*f*t) , 
tl=l/fsl:l/fsl:l; 


signal 
first  ADC 


noisel=sqrt (sigmasq) *randn{l, length(tl) ) ;  %  noise 
ADCsigl=sin(2*pi*f*tl) +noisel;  %  digitized  signal 

ADCsigl=fix(ADCsigl/q) *q;  %  quantized  signal 

t2=l/fs2:l/fs2:l;  %  second  ADC 

noise2=sqrt (sigmasq) *randn(l, length{t2) ) ;  %  noise 
ADCsig2“Sin(2*pi*f*t2) +noise2;  %  digitized  signal 

ADCsig2=fix(ADCsig2/q) *q;  %  quantized  signal 


%f igure (1) 

%subplot (3, 1, 1) ,  plot(t,sig) 
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%title (' Figure  1,  Plot  of  signal') 

%xlabel ( 'Time' ) 

%ylabel ( 'Amplitude' ) 

%subplot (3,1,2) ,plot (tl, ADCsigl) 

%title (' Figure  2.  Plot  of  sampled  signal  (sampling  frequency  1) 
plus  noise') 

%xlabel ( 'time' ) 

%ylabel ( ' magnitude ' ) 

%subplot (3, 1,3), plot (t2, ADCs ig2 ) 

%title(' Figure  3.  Plot  of  sampled  signal  (sampling  frequency  2) 
plus  noise*) 

%xlabel ( ’ time' ) 

%ylabel ( *  magnitude  * ) 

%  Window  operation 
%  Assume  rectangular  window 

winsi2el=fsl; 
winsize2=fs2; 

winsigl=ADCsigl (l:winsizel) ; 
winsig2=ADCsig2 (l:winsize2) ; 

%  DFT  Operation 

DFTsigl=abs (fft (winsigl, winsizel) ) ; 

DFTsig2=abs (fft (winsig2, winsize2) ) ; 

DFTsigla=DFTsigl  (1:  length  (DFTsigl) /2  +1);  %  Taking  half  of  image 
DFTsig2a=DFTsig2  (1:  length  (DFTsig2) /2  +1);  %  Taking  half  of  image 

%figure (2) 

%  Plot  to  locate  position  of  maximum  value 

%  Note  that  due  to  MATLAB  (which  cannot  have  a  zero  index,  the 
actual  location  is  one  less 
%subplot (2, 1, 1) ,  stem(DFTsigla) 

%title  (' Figure  1.  DFT  plot  of  signal  with  sampling  frequency  1*) 
%xlabel  ( *  frequency  bins  ' ) 

%ylabel ( ’magnitude ' ) 

%subplot (2, 1,2) , stem(DFTsig2a) 

%title  (' Figure  2.  DFT  plot  of  signal  with  sampling  frequency  2') 
%xlabel  ( '  frequency  bins ' ) 

%ylabel  ( '  magnitude ' ) 

%  bin  detector 

[i ,yl]  =max (DFTsigla)  ;  %  yl,  y2  are  locations  of  max  values 

[j  ,y2]  =max(DFTsig2a)  ;  %  Note  that  due  to  MATLAB,  the 

%  actual  location  is  one  less. 

al=yl-l; 

a2=y2-l; 

%  SNS  to  Decimal  Algorithm 


%  size  of  window  is  fsl 
%  size  of  window  is  fs2 
%  windowed  sampled  signal  1 
%  windowed  sampled  signal  2 
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%  To  solve  for  f==ai(mod  mi)  (where  "=="  indicates  congruence  and 
%  mi  are  pairwise  relatively  prime) ,  the  Chinese  Remainder  Theorem 
%  states  that  there  is  a  unique  solution  modulo  M=ml*m2 . . . .mr . 

%  A  standard  method  of  solution  is  to  find  integers  hi  such  that 

%  M*bi/mi==l  (mod  mi)  where  i=l,2, _ r  in  which  the  solution  is 

%  f==M*bl*al/ml  +  M*b2*a2/m2  +  ...  +  M*br*ar/mr  (mod  M) 

%  For  a  2  channel  case,  i.e.  i=l,2, 

%  m2*bl  ==  1 (mod  ml) 

%  ml*b2  ==  1 (mod  m2) 

%  f  ==  al (mod  ml) 

%  f  ==  a2 (mod  m2 ) 

%  f  ==  m2*bl*al  +  ml*b2*a2  (mod  ml*m2) 

%  Given  ml (sampling  frequency  1)  and  m2  (sampling  frequency  2),  to 
%  find  bl  and  b2,  the  congruence  equation  is  transformed  to  a 
%  diphantine  equation  and  solved  using  the  Euclidean  algorithm: 

%  m2*bl  -  ml*yl  =  1 

%  ml*b2  -  tn2*y2  =  1 

%  The  above  two  equations  can  be  combined  into 
%  m2*bl  -  ml*b2  =  1 

%  bl  and  b2  are  solved  by  the  function  "Ide.m"  which  is  called  by  % 
"glde.m" . 

% 

%  f  ==  al (mod  ml)  and  f  ==  a2 (mod  m2)  is  solvable  only  if  the 
%  greatest  common  divisor  of  ml  and  m2  divides  (a2  -  al)  . 

%  To  solve  for  f,  r  from  the  diophantine  equation 
%  r*ml+s*m2  =  a2  -  al  must  be  solved. 

%  r  is  obtained  from  "glde.m"  and  f  is  calculated  by  the 
%  equation  f  =  al+r*ml 

idiff=a2-al; 

r  =  glde(fsl,fs2,idiff) ; 
freq=abs (al+r*fsl) ; 

%  Count  the  number  of  correct  results. 

if  freq==f 

count =count + 1 ; 

end 

end 

error  =  1 -count /num; 

%  Write  results  to  file 

xl=fprintf (fpl, '%d  %d  %d  %d  %d  %d\n' ,  f,  fsl,  fs2,  SNRDB,  num,  error); 
plot (SNRDB, error, 'y+') 

title (’Error  Rate  vs.  Signal  to  Noise  Ratio') 
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xlabel ( ' SNR (dB) ' ) 
ylabel { ' Error  Rate  % ' ) 
hold  on 
end 

fclose (fpl) ; 

%  To  calculate  the  dynamic  range 

if  rem(fsl,2) ==0 

DR=fsl/2  +  fs2; 
elseif  rem(fs2,2) ==0 
DR=fs2/2  +  fsl; 

else 

DR=.5* (fsl+fs2) ; 

end 
% 

%  Thesis  Project 
“6 

%  Three  Channel  Receiver 


clear  all; 
close 

%  Initialization 

num= input (* Enter  Number  of  iterations:*);  %  Number  of  iterations 

f= input (* Enter  Input  Frequency:*);  %  Frequency  of  signal 

fsl=input( *  Enter  Sampling  Frequency  1:*);  %  Sampling  frequency  1 

fs2= input  {*  Enter  Sampling  Freqpaency  2:*);  %  Sampling  frequency  2 

fs3=input (’Enter  Sampling  Frequency  3:*);  %  Sampling  frequency  3 

fpl=fopen( * c:\mat lab\bin\thesis\result,dat’ , ’at*) ; 

%  Store  results  in  file  for  later  processing  if  required 

%  Quantization  levels 

%  bit=input ( *  Enter  ADC  resolution: ' ) ; 
bit=14; 

qnlevel=2'^bit“l;  %  No.  of  quantization  levels 

q=2/qnlevel;  %  quantization  size 

for  SNRDB=-30:2:30  %  Set  Signal  to  Noise  Ratio  from  -30 

%  dB  to  30  dB 

count =0; 


%  To  check  whether  fsl  is  even 
%  To  check  whether  fs2  is  even 

%  fsl  and  fs2  are  odd  numbers 


for  i=l:num 

SNR=10^ (SNRDB/10)  ;  %  Convert  to  non-dB  units 

sigmasq=l/2/SNR;  %  Noise  normalization  assuming 

%  signal  power  of  1 


tl=l/fsl:l/fsl:l; 

noisel=sqrt (sigmasq) *randn (1, length (tl) ) ; 

ADCsigl=1000* (sin(2*pi*f*tl) +noisel) ; 

ADCsigl=fix (ADCsigl/q) *q; 

t2=l/fs2:l/fs2:l; 

noise2=sqrt (sigmasq) *randn(l, length (t2) ) ; 

ADCsig2=1000* (sin(2*pi*f *t2) +noise2) ; 

ADCsig2=f ix (ADCsig2/q) *q; 

t3=l/fs3:l/fs3:l; 

noise3=sqrt (sigmasq) *randn(l, length (t3) ) ; 

ADCsig3=1000* (sin(2*pi*f*t3) +noise3) ; 

ADCsig3=f ix (ADCsig3/q) *q; 

%f igure (1) 

^subplot (3,1/1) ,plot (tl,ADCsigl (1 :f si) ) 

%titl€ (* Figure  1.  Plot  of  sampled  signal  (sampling  frequency 
1)  plus  noise') 

%xlabel( 'time' ) 

%ylabel  ( 'magnitude ' ) 

%subplot (3,1,2) ,plot (t2,ADCsig2 (l:fs2) ) 

%t it le (* Figure  2.  Plot  of  sampled  signal  (sampling  frequency  2) 
plus  noise*) 

%xlabel ( 'time* ) 

%ylabel  ( 'magnitude ' ) 

%Subplot (3,1,3) ,plot (t3,ADCsig3 (1 :f s3) ) 

%title (* Figure  3.  Plot  of  sampled  signal  (sampling  frequency  3) 
plus  noise*) 

%xlabel { ' time ' ) 

%ylabel  ( '  magnitude  * ) 


%  first  ADC 
%  noise 

%  digitized  signal 
%  quantized  signal 

%  second  ADC 
%  noise 

%  digitized  signal 
%  quantized  signal 

%  third  ADC 
%  noise 

%  digitized  signal 


%  Window  operation 
%  Assume  rectangular  window 

winsizel=fsl; 
winsize2=fs2; 
winsize3=f s3 ; 

winsigl=ADCsigl (1 rwinsizel) ; 
winsig2=ADCsig2 (1 :winsize2) ; 
winsig3=ADCsig3 (1 :winsize3) ; 

%  DFT  Operation 
DFTsigl=abs (fft (winsigl, wins 


%  size  of  window  is  fsl 
%  size  of  window  is  fs2 
%  size  of  window  is  fs3 
%  windowed  sampled  signal  1 

%  windowed  sampled  signal  2 

%  windowed  sampled  signal  3 

); 


53 


DFTsig2=abs (fft ( wins ig2, wins ize2) ) ; 

DFTsig3=abs (fft (winsig3 , winsize3) ) ; 

DFTsigla=DFTsigl(l: length (DFTsigl)/2  +1);  %  Taking  half  the  image 
DFTsig2a=DFTsig2(l; length (DFTsig2)/2  +1);  %  Taking  half  the  image 
DFTsig3a=DFTsig3(l: length (DFTsig3)/2  +1);  %  Taking  half  the  image 

%f igure (2) 

%Plot  to  locate  position  of  maximum  value 

%Note  that  due  to  MATLAB  (which  cannot  have  a  zero  index,  the 
%actual  location  is  one  less 

%subplot(3,l,l) ,  stem(DFTsigla) 

■s'title  (' Figure  1.  DFT  plot  of  signal  with  sampling  frequency  1') 
%xlabel ( ' frequency  bins ' ) 

%ylabel ( ' magnitude ’ ) 

%subplot(3,l,2) ,stem(DFTsig2a) 

%title(' Figure  2.  DFT  plot  of  signal  with  sampling  frequency  2') 
%xlabel ( ' frequency  bins ' ) 

%ylabel ( ' magnitude ' ) 

%subplot (3,1,3) ,stem(DFTsig3a) 

%t it le (' Figure  3.  DFT  plot  of  signal  with  sampling  frequency  3') 
%xlabel ( ' frequency  bins ' ) 

%ylabel ( 'magnitude ' ) 

%  bin  detector 

[i,yl]=max(DFTsigla) ;  %  yl,  y2  and  y3  are  the  locations  of 

%  maximum  values 

[j .y2] =max(DFTsig2a) ;  %  Note  that  due  to  MATLAB,  the  actual 

%  location  is  one  less. 

[k,y3] =max(DFTsig3a) ; 


al=yl-l 

a2=y2-l 

a3=y3-l 


%  SNS  to  Decimal  Algorithm 

bl=lde(fs2*fs3,fsl) ; 
b2=lde (fsl*fs3, fs2) ; 
b3=lde(fsl*fs2,fs3) ; 

Cl=bl*fs2*fs3; 

C2=b2*fsl*fs3; 

C3=b3*fsl*fs2; 

freqmat= [al*cl+a2*c2+a3*c3;al*cl+a2*c2-a3*c3;  al*cl-a2*c2+a3*c3; 
al*cl-a2*c2-a3*c3;-al*cl+a2*c2+a3*c3;-al*cl+a2*c2-a3*c3; 
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-al*cl-a2*c2+a3*c3; -al*cl“a2*c2-a3*c3]  ; 

f reqmat=rem ( f reqmat , f sl*f s2*f s3 ) ; 

for  i=l:8 

if  (freqmat (i) <0) 

freqmat (i) =freqmat (i) +f sl*fs2*fs3 ; 

end 

end 

freq=min(abs  (freqmat) )  ; 

%  Count  the  number  of  correct  results. 

if  freq^=f 

count =count+l ; 

end 

end 

error  =  1 -count /num; 

%  Write  results  to  file 

xl=fprintf (fpl, * %d  %d  %d  %d  %d  %d  %d\n* ,  f,  fsl,  fs2,  fs3/  SNRDB,  num, 
error) ; 

plot (SNRDB, error, *y+*) 

hold  on 

end 

fclose (fpl) ; 

%  To  calculate  the  dynamic  range 

if  rem(fsl,2) ==0  %  To  check  whether  fsl  is  even 

x=[fsl/2  +  fs2*fs3;  fsl*fs2/2  +  fs3;  fsl*fs3  +  fs2] ; 
elseif  rem(fs2,2) “=0  %  To  check  whether  fs2  is  even 

x=[fs2/2  +  fsl*fs3;  fsl*fs2/2  +  fs3;  fs2*fs3  +  fsl] ; 
elseif  rem(fs3, 2) ==0  %  To  check  whether  fs3  is  even 

x=[fs3/2  4*  fs2*fsl;  fs3*fs2/2  +  fsl;  fsl*fs3  +  fs2]  ; 
else  %  fsl,fs2  and  fs3  are  odd 

x=l/2* [f si  +  fs2*fs3;  fs2  +  fsl*fs3;  fs3  +  fs2*fsl]  ; 

end 

DR=min  (x)  ; 


55 


%  This  function  solves  the  general  linear  diophantine  equation 
%  in2*bl  -  ml*b2  =  k  and  returns  the  value  bl 

function  a=glde (ml,m2,k) 

%  Calls  function  "Ide”  to  calculate  bl,  b2  and  na 
[bl , b2 , na] =lde (ml , m2 ) ; 

%  To  check  whether  the  equation  is  solvable. 

%  na  must  be  a  factor  of  k  for  the  equation  to  be  solvable. 

mult=k/na; 

if  {k-mult*na) ==0  %  Equation  is  solvable 

bl=bl*mult;  %  These  new  values  solve  the  diophantine  equation 
b2=b2*mult ; 

mtest=bl;  %  To  check  whether  bl  and  b2  are  the  least  values 

mdl=ml/na;  %  that  satisfies  the  diophantine  equation 

md2=m2/na; 

mx=bl; 

mx=mx+md2 ; 

while  (abs (mx) -abs (bl) ) <0 
bl=mx; 
b2=b2”mdl; 
mx=tnx+md2 ; 

end 

if  (mtest-bl) ==0 
mx=bl; 
mx=mx-md2; 

while  (abs (mx) -abs (bl) ) <0 
bl=mx; 
b2=b2+mdl; 
mx=Tnx-md2  ; 

end 

end 

end 
a=bl  ; 
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%  This  function  solves  the  linear  diophantine  equation 
%  ml*bl  +  m2*b2  =  na  where  ml  and  m2  are  the  sampling  frequencies 
%  and  na  is  the  greatest  common  divisor 
%  and  returns  the  value  bl,  b2  and  na 
% 

%  ml  and  m2  are  assumed  positive 

function  [bl , b2 , na] =lde (ml , m2 ) 

%  Initialize  bol,  bo2,  bl  and  b2 

bol=l; 

bo2=0; 

bl=0; 

b2=l; 

%  Place  ml  and  m2  in  ma (dividend)  and  na  (divisor)  respectively 

ma=ml ; 
na=m2 ; 

%  Calculate  quotient  and  remainder 

iquot-f ix (ma/na) ; 
irem=ma-na*iquot ; 

%  If  remainder  is  not  zero,  reset  dividend  and  divisor 

while  irem>0 

bo  3  =bo 1 - iquo t  *b 1 ; 
bo4  =bo2 - iquot *b2 ; 
bol=bl; 
bo2=b2; 
bl“bo3 ; 
b2=bo4; 
ma=na; 
na=irem; 

ic[uot=fix  (ma/na)  ; 
irem=ma-na*iq[uot  ; 

end 


%  calculate  new  coefficients  of  ml  and  m2 
%  redefine  bol,  bo2,  bl  and  b2 

%  redefine  dividend  and  divisor 
%  reapply  Euclidean  algorithm 
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TMS320C54X  DSKplus 

AJeebraic  Assembler 


Ke\'  Features 


Code  Explorer  Debugger 
System  Requirements 
How  to  Install 


The  C54x  DSKplus  is  a  low-cost  design  tool  that  gives  designers  a 
working  knowledge  of  DSP  code.  From  this  foundation,  designers  can 
begin  buflding  complete  'C54x  DSP-based  systems.  Priced  at  US  $149, 
the  •C54x  DSKplus  (part  no.  TMDS32000L0)  is  available  from  Tl 
airthorized  distributors . 


Bevond  the  DSKplus 
*C54x  Software  Sunoort  Files 


The  C54x  DSKplus  builds  on  TTs  industry-leading  line  of  low  cost, 
e^-to-use  DSP  Starter  Kit  PSK)  development  boards.  The 
high-performance  board  features  the  TMS320C542  16-bit  fixed-point 
DSP.  Capable  of  performing  40  million  instructions  per  second 
(MIPS),  the  'C542  makes  the  'C54x  DSKplus  the  most  powerful  DSK 
development  board  on  the  market. 

Other  TMS320  DSKs  include  the  'C2xDSK.  the  *C5x  DSK .  and  the 
floating-point  ‘C3x  DSK . 


Key  Features 

The  'C54x  DSKplus  includes: 

•  40  MIPS  TMS320C542-b3sed  board 

•  TLC320AC01  Analog  Interfece  Circuit  (AIC) 

•  'C54x  DSKplus  assembler,  loader.  Code  Explorer  debugger,  and  sample  programs  (3.5"  disks) 

•  TMS320C54X  CPU  and  Peripherals  Reference  Guide 

•  TMS320C54X  Algebraic  Assembler  Instruction  Set 

•  TMS320C54X  Datasheet 

•  TMS320C54X  DSKplus  User’s  Guide 

•  TLC320AC01  Datasheet 

•  PC  coimector  cable  and  universal  power  supply  included 

•  US‘5149  discount  coupon  toward  the  purchase  of  the  'C54x  EVM 
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DSKplus 

Features  I 

Benefits  j 

TMS320C542  DSP  ^0 

MIPS,  16-bit) 

High-performance,  very  efScient  architecture  requires  fewer  | 

MIPS  than  conipeting  DSPs  to  implement  most  algorithms.  ? 

Code  Explorer  debugger  | 

interface 

An  easy-to-use,  true  Windows-based  interface.  Supports  | 

symbolic  debugging,  breakpoints,  graphical  animation,  variable  [ 
watch  windows,  JBle  FO,  algebraic/mnemonic  disassembfy,  j 

on-line  help.  1 

Symbolic  debugging  (Code  \ 
Explorer) 

. 1 

1 

Enables  easy  programmability  by  using  labels  for  referencing  \ 

constants,  variables,  matrices  name.  1 

.  Algebraic  assembler 

1 

Bypasses  learning  new  DSP  mnemonic  instruction  set  specifics.  : 
Makes  coding  easier  and  more  straight-forward.  Easy  one-step  j 

assembly  and  linking  process.  1 

j 

Demo  programs/ 

Application  code 

j  1 

i 

:  Helps  users  get  up-to-speed  quickly  t 

i 

.  J 

TLC320AC01  Analog 
Interface  Chip 

Low  power  dissipation,  14-bit  linear  resolution,  programmable  ; 

sampling  rates,  anti-aliasing  filter,  and  input  gain;  selectable  j 

auxQiaiy  input;  data  read-back  1 

> 

j 

Socketed  Programmable 
Array  Logic  (PAL) 

Allows  e.^qjerienced  designers  to  reprogram  the  PAL  and  change  1 
the  w^  the  host  port  interface  works  on  the  C54x  DSKplus.  j 

Universal  power  supply  & 
cable  included  : 

i 

Allows  for  immediate  use  out  of  the  box;  ideal  for  powering  \ 

daughter  cards;  filtered  and  regulated  -  finis  no  need  for  | 

on-board  voltage  regulation.  1 

4 

’C54x  Algebraic  Assembler 

The  CS4x  DSKplus  includes  the  algebraic  assembler  that  speeds  the  initial  code  development  process.  The 
algebraic  assembler  does  not  require  new  users  to  learn  a  new  DSP  mnemonic  instruction  set,  maVing 
coding  easier  and  more  direct  The  assembler  also  utilizes  a  one-step  assembly  and  linking  process  to 
simplify  code  ddwgging.  The  software  accomplishes  this  by  using  special  directives  to  assemble  code  at  an 
absolute  address. 

Some  extremefy  useful  features  include: 
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•  In-line  Assembly  expression  analysis  ^ows  the  assembler  to  work  when  defining  complex  variables 
orbit  locations. 

Symbolic  Debugg^g  allows  the  user  to  reference  variables  by  name  instead  of  the  physical  address. 

•  Assembling  conditional  blocks  of  assembly  code  using  .if'.else  if'.end  if  directives.  This  is  espedalty 
helpful  when  you  want  to  conditionally  assemble  code  via  a  command-line  argument  of  internal 
assembly  variable. 

•  Support  of  .sect,  .bss,  .usect,  .text,  and  .data  sections. 


Code  Explorer  Debugger 


The  •C54x  DSKplus  debugger  was  developed  by  GO  DSP  Corporation  in  an  effort  to  provide  the  first  true 
Windows-based  debugger  for  a  DSK.  The  Code  Explorer  debugger  siq5X)rts  debugging,  a  new  feature 
available  only  on  the  DSKplus  that  allows  the  user  to  specify  labels  for  referencing  constants,  variables, 
and  marticies  ly  name.  Also,  the  debugger  desktop  environment  is  fully  configurable  and  loaded  iqwn  * 
entry  into  the  debugger.  This  means  that  optional  colors,  fonts,  and  window  sizes  can  be  within 

the  dd>ugger  and  saved  upon  exiting. 


Some  additional  features  of  the  debugger  include  capability  of  connecting  files  as  I/O,  graphical  animatinn 
and  ^ta  memory  viewing.  The  file  I/O  capability  enables  users  to  connect  files  as  inputs  or  outputs  to  aity 
location  within  your  application  code.  Therefore  you  can  simulate  different  input  setyiences  and  data 
streams  without  having  to  pitysically  generate  them. 

Graphical  animation  allows  you  to  view  data  in  a  graphical  format,  either  with  time  HnmaiTi  or  frequency 
domain  and  in  a  variety  of  variable  sizes  (i.e.  8-bit  signed  char,  8-bit  unsigned  char,  16-bit,  32.bit,  etc). 

Disassembly  Window 


The  disassembly  window  displays  the  DSP  code  in 
algebraic  instructions.  The  variable  names  and 
subroutines  (symbols)  are  shown  in  blue.  The 
physical  DSP  address  is  the  first  column  and  the 
machine  code  for  the  instructions  are  in  column  2. 
The  yellow  bar  indicates  the  location  where  the 
DSP  program  counter  (PC)  points. 

The  disassembly  window  properties  can  be 
accessed  by  placing  the  cursor  in  the  disassembly 
window  and  right-clicking  and  then  choosing 
properties.  The  disassembly  window  can  display 
code  in  algebraic  or  mnemonic  formats  with  direct 
and  immediate  addressing  values  shown  in  hex, 
decimal  and  even  binary. 

Data  Memory  Window 
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The  data  memory  window  can  be  modified  or 
replicated  as  needed.  By  placing  the  cursor  inside 
the  data  memory  window  and  right-clicking  and 
then  choosing  properties,  the  user  can  change  the 
title  of  the  window,  starting  address  and  even  data 
organization  in  the  window.  Valid  display  formats 
include  8-bit  signed/unsigned  char, 
signed/unsigned  long,  floats,  and  others.  The  page 
field  can  q)ecify  either  Data  or  Program  memory 
spaces. 

’C54x  CPU  and  Peripheral 


Registers 


The  two  register  windows  in  the  'C54x  Code 
Ejq5lorer  debugger  are  the  CPU  and  Peripheral 
Registers.  The  'C54x  CPU  Registers  is  the 
collection  of  registers  which  control  the  operation 
of  the  DSP  CPU.  The  program  cormter,  status 
register,  and  configuration  re^sters  are  contained 
within  this  window.  Notice  that  bit  values  within 
the  register  are  brought  out  separately  to  make 
modification  and  monitoring  easier. 

The  second  window  is  the  Peripherals  window. 
This  window  includes  the  registers  for  configuring 
the  DSP  peripherals  like  the  serial  ports  and 
timers.  Modifications  to  this  register  can  be  done 
by  clicking  on  the  register  in  the  Peripheral 
Registers  window. 


Graphical  Windows 


Graphical  windows  are  extremely  useful 
when  trying  to  view  a  value  of  a  register, 
variable,  or  buffer.  The  graphic  window 
allows  the  user  to  animinate  any  value  in 
either  data  or  program  DSP  memory. 
This  is  accomplished  by  placing  a 
breakpoint  anywhere  in  the  application 
code  and  pressing  the  Animation  button. 
Each  time  the  DSP  reaches  the 
brealqx)int  the  graphical  windows  are 
updated  and  refieshed. 

The  options  window  contains  the 
graphics  setup  for  the  window.  For 
example,  the  title  can  be  changed  to 
reflect  the  data  being  animated,  the 
display  buffer  length  can  be  changed,  or 
the  data  read  firom  the  DSP  can  either  be 


a  qngift  value  from  a  list  (buffer)  of  values  in  either  data  or  program  memory.  Also,  the  sampling  rate  can 
be  modified  for  correct  diqjlaying  of  the  frequency  data  (FFT).  The  diq)lay  can  be  viewed  using  8-bit 
signed/unsigned  chars,  infs,  long  floats,  and  even  a  log  can  be  performed  on  the  displayed  data. 


Setting  Breakpoints 
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A  brealqxjint  can  be  selected  by  either 
double  clicking  on  a  line  in  the 
Disassembly  window  or  by  Selecting  the 
DEBUG-BREAKPOINTS  in  the  PuU 
Down  Menu.  The  Pull  Down  Menu  will 
pronqjt  you  with  a  menu  listing  all  the 
available  symbols  in  the  Symbols  box. 
You  can  either  select  a  brealqx>int  from 
the  list  of  Symbols  or  by  entering  an 
address  in  the  Address  field. 


The  Brea]qx)int  dialog  box  contains  the 
following  fields:  Address,  Symbols  and 
Breakpoints.  If  the  address  of  the  desired 
breakpoint  is  known,  simpfy  enter  the 
value  in  the  Address  field.  The  Symbol 

field  contains  the  list  of  afi  the  symbols  in  the  program.  If  the  location  address  of  the  breakpoint  is  labeled. 
Simply  type  the  label  name  and  press  add 

Setting  Probe  Points 

Probe  points  allow  the  update  of  a 
particular  window  or  the  reading/writing 
of  samples  from  a  file  to  occur  at  a 
specific  point  in  an  algorithm.  This 
effectively  **connects  a  signal  probe"  to 
that  point  in  the  algorithm. 


When  a  graph  window  otgect  is  created, 
it  assumes  that  it  is  to  be  updated  at  every 
hrealqx)int.  However,  this  attribute  can  be 
changed  and  the  window  can  be  updated 
only  when  the  program  reaches  the 
connected  probe  point.  After  the  probe 
point  is  hit,  and  the  window  is  updated, 
execution  of  the  program  is  continued 
This  optimizes  the  displ^  of  the  graph 
window  and  also  allows  you  to  keep  a  history  of  the  signal  even  when  the  data  on  the  DSP  is  not  valid 

With  the  combination  of  Code  Explorer's  File  I/O  cap^ilities,  probe  points  can  also  used  to  connect 
streams  of  data  to  a  particular  point  in  the  DSP  Code.  When  the  probe  point  is  reached  in  the  algorithm, 
data  is  streamed  from  a  specific  memory  area  to  file,  or  from  tire  file  to  memory. 

Using  File  I/O 


I 
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Code  E^lorer  allows  the  user  to  stream 
data  onto  (or  from)  the  target  from  a  PC 
file.  This  ^ows  the  user  to  simulate  code 
tising  known  sample  values.  Note  that 
this  file  I/O  feature  is  not  intended  to 
satisfy  real-time  constraints.  The  File 
Input/Output  feature  uses  probe  points. 

When  the  execution  of  the  program 
reaches  a  probe  point,  the  connected 
olyect,  whether  it  is  a  file,  graph  or 
memoiy  window,  is  updated.  Chice  the 
connected  object  is  igxlated,  execution 
continues.  Using  this  concept,  if  a  probe 
point  is  set  at  a  specific  point  in  the  axle 
and  then  connected  to  a  file,  file  I/O  functionalities  can  be  implemented 


System  Requirements 

•  A386,486,  or  Pentium  PC  with  a  3.5"  disk  drive 

•  4-bit  parallel  and/or  843it  bidirectional  parallel  ports. 

•  A  minimum  of  4Mhytes  of  memory 

•  Color  VGA  monitor 

•  Windows  3.1  or  Windows  95 

•  ASCn  editor 


How  to  Install 

When  connecting  the  DSKplus  to  your  PC,  it  is  hi^y  recommended  you  turn 
off  your  PC’s  power  to  make  the  coimections  below: 


1.  Connect  the  DB25  cable  (female)  to  the  PC's  Parallel  port  (male). 

2.  Cormect  the  DB25  cable  (male)  to  the  DSKplus  board  (female). 

3.  Cormect  the  power  cord  (NEMA  cable)  to  the  5  volt  power  suppfy. 

4.  Connect  the  5-pin  DIN-to-5.5mm  adapter  to  the  power  supply’s  5-pin 
DIN  coimector. 

5.  Plug  the  power  suj^ly  power  cord  to  the  wall  outlet 

6.  Plug  the  5.5mm  coimector  into  the  power  jack  of  the  DSKplus  board 

At  this  point  the  green  power  LED  is  illuminated  and  power  is  siq^lied  to  the  ’C54x  DSKplus  board  If  the 
Green  LED  is  not  iUuminated,  check  the  coimections  on  the  power  si^fy  and  power  cord 

Installing  the  software 
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TTie  DSKplus  kit  includes  two  3.5"  floRjies  l^led  Disk  #1  and  Disk  #2.  To 
install  the  software  correctly,  please  follow  the  steps  below; 

1.  Insert  Disk  #1  into  the  3.5"  drive. 

2.  From  the  start  menu  (Windows95)  or  the  Files  menu  (Windows  3.1) 
select  the  Run.,  oj^oa  Type  a:\setup.exe 

3.  The  installation  script  will  af^)ear.  You  will  be  asked  to  select  a 
destination  directory.  By  default  it  will  select  the  DSKplus  directory. 

Enter  the  directory  name  if  you  would  like  to  q^ecify  a  rfiffprpnt 
directory. 

4.  When  prompted,  insert  Disk  #2  into  the  3.5"  floppy  drive. 

5.  When  installation  has  oorrqrleted,  the  installation  will  irrform  you  that  the  installation  was 
successful.  At  this  point  a  Code  Erqjlorer  Grorq)  will  a^jear. 


Starting  the  Debugger 

To  start  the  ^bugger,  click  on  the  icon  located  in  the  Code  Ejqjlorer  Group  or  deddop.  The  Code  E^gjloier 
background  and  winctows  will  ^^sear  with  the  Setup  Box  shown  active. 


Select  the  port  which  is  cormccted  to  the  DSKplus  board.  If  for  some  reason 
the  port  is  not  listed,  the  port  address  can  be  modihedl^  typing  in  the  address 
int  the  text  box. 

As  a  result  of  selecting  the  correct  port  and  proper  hardware  coimections,  the 
debugger  will  fill  its  windows  with  data  and  the  DSKplus  is  now  fimrtinning 
If  for  some  reason  the  ddnigger  responds  with  the  error  "Can't  initialiTP 
Target  DSP",  follow  the  directions  in  the  error  box. 


Troubleshooting 


1.  Is  the  power  on?  Be  sure  green  LED  is  illuminated.  If  not,  a  loose  power  cable  is  hampering  your 
setip. 

2.  Is  the  panllel  port  cable  connection  secure?  In  many  new  DSKplus  boards  and  parallel  port 
substantial  pressure  maity  be  needed  to  connect  the  cables.  Connect  the  cable  to  the  DSKjdus  board 
by  placing  the  thumb  behind  the  DB-25  connector.  Take  the  caWe  connector  chassis  and  place 
between  the  index  and  middle  fingers.  Align  the  connectors  and  press  &e  fingers  together. 

3.  The  port  selected  is  not  being  "Ctgjtured"  by  Windows  95.  Capturing  is  used  by  Windows  95  to 
allow  DOS  programs  access  to  printers.  The  port  can  be  released  by  going  into  the  control  panel 
and  selecting  the  printers  icoa  Highlight  any  printer  and  go  to  the  File  pulldown  on  the  command 
bar.  Select  properties  and  then  the  Details  tab.  The  Details  tab  includes  a  button  nanwvt  End 
Capture...  Click  on  this  bitton  and  select  the  LPT  port  where  the  DSKplus  board  is  connected.  If 
the  LPT  is  not  listed,  then  the  port  is  not  captured  (select  cancel)  and  proceed  to  number  4. 

4.  The  port  sdected  is  configured  as  an  EPP  or  ECP  port.  The  DS^lus  board  si^qxjrts  4-bit 
unidirectional  and  8-bit  bidirectional  parallel  ports.  The  DSKplus  does  not  support  EPP  and  ECP 
ports.  To  check  the  port  configuration,  «dt  out  and  rdxx)t  your  system.  At  the  point  where  the 
BIOS  Setup  routine  can  be  selected,  press  the  keyboard  sequence  to  enter  the  BIOS  (usually 
CTRL+ALT+ESC).  Confirm  that  the  parallel  port  is  setup  as  '8-bit’,  'bidirectional’  or  'standard.' 
Specifically,  not  an  EPP  or  ECP  port  If  problms  persist,  run  the  included  selftest  program. 


Beyond  the  ’C54x  DSKplus 

With  higher  performance  than  any  other  DSK  available  to(%,  the  'C54x  DSKplus  offers  a  rich 
development  enviromnent  for  benchmarking  and  evaluating  code  in  real-time.  The  ’C54x  DSKplus  is 
designed  as  an  eaq^-to-use  entry  into  the  worid  of  high-performance  fiixed-point  DSPs. 
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However,  as  your  design  e^qjeiience  grows,  you  require  addition^  functionality  and  expanded 
c^abilMes.  To  meet  these  needs,  TI  offers  a  conqjrehensive  line  of  evelopment  tools  for  the  TMS320  DSPs 
that  Siq^rt  the  design  process  from  system  concept  to  production. 

DtTier  *C54x  PeveloOTient  Tools 


©  Copvrishi  1997  Texas  Instruments  Incorporated.  All  ri^its  reserved 
Trademarks.  Important  Notice! 
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APPENDIX  C 


C  AND  ASSEMBLY  LANGUAGE  CODE  FOR  DSP  HARDWARE 


**********■***************************★*******-*****★****•**************** 
File:  FirstApl.ASM 

When  sampling  frequency  is  changed,  need  to  change 

a.  buffer  size  here 

b.  A  and  B  registers  in  ACOlinil.asm 

c.  sampling  frequency  in  dftsort.c 

d.  buffer  size  in  hostappl.cpp 

★*******************•****•★*****************'***********’***************★** 
.width  80 
.length  55 

.title  "FirstApp  program” 

.mmregs 

.setsect  ".text",  0x500,0 
.setsect  "vectors",  0x180,0 


VECTORS 


.sect  "vectors" 

. copy  "c : \dskplus\inits\vectors . asm" 
.text 


start : 

call  ACOIINIT 
pmst  =  #01a0h 
sp  =  #0ffah 
ar2  =  #1200h 
*ar2+  =  data(#0bh) 
imr  =  #280h 
intm  =  0 

wait  nop 

goto  wait 


;  set  up  iptr 
;  init  stack  pointer. 

;  pointer  to  receive  buffer  at  120 Oh. 
;  store  to  rev  buffer 

;  ready  to  rev  int ' s 


XINT: 


Receive  Interrupt  Routine 


b  -  trev  ; 

b  =  #0FFFCh  &  b 
*ar2+  =  data(#0bh)  ; 

tdxr  =  b  ; 

TC  =  (@ar2  ==  #1471h)  ; 

if  (TC)  goto  restrt 
return  enable 


load  acc  b  with  input 

store  to  rev  buffer 
transmit  the  data, 
change  here  if  fs  changes 
stop  if  rev  buffer  is  at  1471h 


restrt 
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ar2  =  #1200h 
hpic  =  #Oah 


;  set  intm  bit  ...no  int* *s 
;  flag  host  task  completed 


return  enable 


end  I SR 


. copy  ”c : \dskplus\f irstapp\ac01inil . asm" 
.end 


*********************************************************************** 
File:  FirstApp.ASM  ->  First  Application  program  for  the  *C54x  DSKplus 

a.  buffer  size  here 

b.  A  and  B  registers  in  AC01ini2.asm 

c.  sampling  frequency  in  dftsort.c 

d.  buffer  size  in  hostapp2.cpp 


.width  80 
.length  55 

.title  "FirstApp  program" 

. mmregs 

.setsect  ".text",  0x500,0 
.set sect  "vectors",  0x180,0 


VECTORS 


.sect  "vectors" 

.copy  "c:\dskplus\inits\vectors.asm" 

.text 

start : 

call  ACOlINIT 

pmst  =  #01a0h  ;  set  up  iptr 

sp  =  #0ffah  ;  init  stack  pointer. 

ar2  =  #l200h  ;  pointer  to  receive  buffer  at  1200h. 

*ar2+  =  data (#0bh)  ;  store  to  rev  buffer 

imr  #280h 

intm  =0  ;  ready  to  rev  int ’ s 


wait  nop 

goto  wait  . 


Receive  Interrupt  Routine 


XINT: 

b  =  trev 
b  =  #0FFFCh  &  b 
*ar2+  =  data(#0bh) 
tdxr  =  b 

TC  =  {@ar2  ==  #01400h) 
if  (TC)  goto  restrt 
return  enable 


;  load  acc  b  with  input 

;  Store  to  rev  buffer 
;  transmit  the  data. 

;  change  here  if  fs  change 
;  stop  if  rev  buffer  is  at  1400h 


restrt 


;  set  intm  bit  ...no  int's 
;  flag  host  task  completed 


ar2  =  #120 Oh 
hpic  =  #0ah 
return_enable 

end  ISR 

•copy  "C:\dskplus\firstapp\ac0lini2.asm" 

.end 


************** 


****************** 


File:  ACOIINII.ASM  ->  ACOl  Initialization  Routine 

*********************************************************************** 

.width  80 
.length  55 

.title  "ACOl  Initialization  Program" 

.mmregs 

************************************************************************ 
**Certain  ACOl  registers  can  be  initialized  using  a  conditional  assenibly 
constant.  By  setting  the  constant  REGISTER  to  the  appropriate  value, 
thrLsemblL  will  either  include  initialization  for  certain  registers 
or  ignore  register  initialization. 

The  constant  REGISTER  should  be  set  to  include  the  following  ACOl 
register : 


* 

REGISTER 

(binary)  = 

* 

0000 

0000 

0000 

0001 

* 

0000 

0000 

0000 

0010 

ic 

0000 

0000 

0000 

0100 

* 

0000 

0000 

0000 

1000 

ie 

* 

0000 

0000 

0001 

0000 

* 

ie 

0000 

0000 

0010 

0000 

★ 

* 

0000 

0000 

0100 

0000 

ic 

0000 

0000 

1000 

0000 

->  initialize  Register  1 
->  initialize  Register  2 
->  initialize  Register  3 
->  initialize  Register  4 

->  initialize  Register  5 

->  initialize  Register  6 

->  initialize  Register  7 
->  initialize  Register  8 


(A  Register) 

(B  Register) 

(A'  Register) 
(Amplifier  Gain- 
Select) 

(Analog 

Conf igurat ion) 
(Digital 
Conf igurat ion) 
(Frame-Sync  Delay) 
(Fram-Sync  number) 


Any  combination  of  registers  can  be  initialized  by  adding  the  binary 
number  to  the  REGISTER  constant.  For  example  to  initalize  Registers 
and  5,  REGISTER  =  18h.  Upon  assembly,  only  code  for  register  4  &  5 
initialization  is  included  in  the  ACOIINIT  module.  When  called  the 
module  will  load  REG4  and  REGS  values  into  internal  ACOl  registers. 


*  Register  4  is  always  loaded  to  get  a  6db  input  gain.  This  sets  full- 

*  sSle  to  3v(p-p  input)  due  to  the  single-ended  ACOl  configuration. 


REGISTER  .set 


REGl 

REG2 

REG3 

REG4 

REGS 

REG6 


.set 

.set 

.set 

.set 

.set 

.set 


Obh 

If  eh 

21fh 

300h 

40dh 

SOlh 

600h 


Powerup  default  values ; 

*  112h 

*  212h 
300h 

*  40Sh 

SOlh 

600h 
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REG7  .set  700h 

REGS  .set  SOlh 


700h 

SOlh 


AGO  UNIT: 

xf  =  0 
intm  =  1 
ter  =  #10h 
imr  =  #280h 
tspc  =  #000Sh 
tdxr  =  #0h 
tspc  =  #00cSh 
xf  =  1 


reset  acOl 

disable  all  int  service  routines 
stop  timer 

wakeup  from  idle  when  TDM  Xmt  int 
stop  TDM  serial  port 
send  0  as  first  xmit  word 
reset  and  start  TDM  serial  port 
release  acOl  from  reset 


Register  init ' s 


.eval  REGISTER  &  Ih,  SELECT  ;  if  REGl  then  include  this  source 
.if  SELECT  =  Ih 

a  =  #REG1  ;  load  Acc  A  with  REGl  value 

call  REQ2  ;  Call  REQ2  subroutine 

.endif 

.eval  REGISTER  &  2h,  SELECT  ;  if  REG2  then  include  this  source 
.if  SELECT  =  2h 
a  =  #REG2 
call  REQ2 


.endif 

.eval  REGISTER  &  4h,  SELECT  ;  if  REGS  then  include  this  source 

.if  SELECT  =  4h 

a  =  #REG3 

call  REQ2 

.endif 

.eval  REGISTER  &  Sh,  SELECT  ;  if  REG4  then  include  this  source 
.if  SELECT  =  Sh 
a  =  #REG4 
call  REQ2 
. endif 

.eval  REGISTER  &  lOh,  SELECT  ;  if  REGS  then  include  this  source 

.if  SELECT  =  lOh 

a  =  #REG5 

call  REQ2 

.endif 

.eval  REGISTER  &  2 Oh,  SELECT  ;  if  REG6  then  include  this  source 

.if  SELECT  =  20h 

a  =  #REG6 

call  REQ2 

.endif 
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.eval  REGISTER  &  4 Oh,  SELECT  ;  if  REG7  then  include  this  source 
•if  SELECT  =  40h 
a  =  #REG7 
call  REQ2 
. endif 


.eval  REGISTER  &  80h,  SELECT  ;  if  REGS  then  include  this  source 

•if  SELECT  =  80h 

a  =  #REG8 

call  REQ2 

. endif 

return 


REQ2 

ifr  =  #080h  ;  clear  flag  from  IFR 

tdxr  =  #03h  ;  request  secondary  when  ACOl  starts 


idled) 
tdxr  =  a 
ifr  =  #080h 


;  wait  for  primary  to  xmit 
;  send  register  value  to  serial  port 
;  clear  flag  from  IFR 


idled) 
tdxr  =  #0h 
ifr  =  #080h 
idled) 
return 
.end 


;  wait  for  secondary  to  xmit 
;  send  neutral  state  in  case  last  init 
;  clear  flag  from  IFR 
;  wait  for  neutral  state  to  xmit 
;  return  from  subroutine 
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-*********************************************************************** 

/ 

7 

;  File:  AC01INI2 .ASM  ->  ACOl  Initialization  Routine 


.*********************************************************************** 

/ 


.width  80 
.length  55 

.title  ”AC01  Initialization  Program" 
.mmregs 


************************************************************************ 


*  Certain  ACOl  registers  can  be  initialized  using  a  conditional  assembly 

*  constant.  By  setting  the  constant  REGISTER  to  the  appropriate  value, 

*  the  assembler  will  either  include  initialization  for  certain  registers 

*  or  ignore  register  initialization. 

* 


*  The  constant  REGISTER  should  be  set  to  include  the  following  ACOl 

*  register: 

* 


* 

ie 

REGISTER 

(binary) 

= 

* 

0000 

0000 

0000 

0001 

* 

0000 

0000 

0000 

0010 

* 

0000 

0000 

0000 

0100 

* 

ie 

0000 

0000 

0000 

1000 

* 

it 

0000 

0000 

0001 

0000 

* 

ie 

0000 

0000 

0010 

0000 

* 

0000 

0000 

0100 

0000 

* 

0000 

0000 

1000 

0000 

“>  initialize  Register  1 
-">  initialize  Register  2 
->  initialize  Register  3 
->  initialize  Register  4 

->  initialize  Register  5 

->  initialize  Register  6 

->  initialize  Register  7 
->  initialize  Register  8 


(A  Register) 

(B  Register) 

(A’  Register) 
(Amplifier  Gain- 
Select) 

(Analog 

Conf igurat ion ) 
(Digital 
Conf igurat ion) 
(Frame-Sync  Delay) 
(Fram-Sync  number) 


*  Any  combination  of  registers  can  be  initialized  by  adding  the  binary 

*  number  to  the  REGISTER  constant.  For  example  to  initalize  Registers  4 

*  and  5,  REGISTER  =  18h,  Upon  assembly,  only  code  for  register  4  &  5 

*  initialization  is  included  in  the  ACOIINIT  module.  When  called  the 

*  module  will  load  REG4  and  REG5  values  into  internal  ACOl  registers. 


* 


* 

*  Register  4  is  always  loaded  to  get  a  6db  input  gain.  This  sets  full- 


*  scale 

ie 

to  3v(p 

-p  input) 

due  to  the  single 

-ended  ACOl  configuration 

ie 

REGISTER 

.set 

Obh 

;  Powerup  default 

values : 

REGl 

.set 

ifeh 

.  * 

/ 

112h 

REG2 

.set 

23ch 

-  * 

/ 

212h 

REG3 

.set 

300h 

7 

300h 

REG4 

.set 

40dh 

.  * 

/ 

405h 

REG5 

.set 

501h 

7 

501h 

REG6 

.set 

600h 

r 

600h 
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REG7 

REGS 


.set  700h 
.set  SOlh 


700h 

SOlh 


AGO  UNIT: 

xf  =  0 
intm  =  1 
ter  =  #10h 
imr  =  #2S0h 
tspc  =  #OOOSh 
tdxr  =  #0h 
tspc  =  #00cSh 
xf  =  1 


reset  acOl 

disable  all  int  service  routines 
stop  timer 

wakeup  from  idle  when  TDM  Xmt  int 
stop  TDM  serial  port 
send  0  as  first  xmit  word 
reset  and  start  TDM  serial  port 
release  acOl  from  reset 


-  Register  init's  - 

.eval  REGISTER  &  Ih,  SELECT  ;  if  REGl  then  include  this  source 
.if  SELECT  =  Ih 

a  =  #REG1  ;  load  Acc  A  with  REGl  value 

call  REQ2  ;  Call  REQ2  subroutine 

.endif 

.eval  REGISTER  &  2h,  SELECT  ;  if  REG2  then  include  this  source 
.if  SELECT  =  2h 
a  =  #REG2 
call  REQ2 

.endif 

.eval  REGISTER  &  4h,  SELECT  ;  if  REG3  then  include  this  source 

.if  SELECT  =  4h 

a  =  #REG3 

call  REQ2 

.endif 

.eval  REGISTER  &  8h,  SELECT  ;  if  REG4  then  include  this  source 

.if  SELECT  =  8h 

a  =  #REG4 

call  REQ2 

.endif 

.eval  REGISTER  &  lOh,  SELECT  ;  if  REGS  then  include  this  source 

.if  SELECT  =  lOh 

a  =  #REG5 

call  REQ2 

.endif 

.eval  REGISTER  &  2 Oh,  SELECT  ;  if  REG6  then  include  this  source 
.if  SELECT  =  20h 
a  =  #REG6 
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call  REQ2 
.endif 

.eval  REGISTER  &  40h,  SELECT  ;  if  REG7  then  include  this  source 

.if  SELECT  =  40h 

a  =  #REG7 

call  REQ2 

.endif 


.eval  REGISTER  &  8 Oh,  SELECT  ;  if  REGS  then  include  this  source 

.if  SELECT  =  80h 

a  =  #REG8 

call  REQ2 

.endif 

return 


REQ2 

ifr  =  #080h  ;  clear  flag  from  IFR 

tdxr  =  #03h  ;  request  secondary  when  ACOl  starts 


idled) 
tdxr  =  a 
ifr  =  #080h 


;  wait  for  primary  to  xmit 
;  send  register  value  to  serial  port 
;  clear  flag  from  IFR 


idled) 
tdxr  =  #0h 
ifr  =  #080h 
idled) 
return 
.end 


wait  for  secondary  to  xmit 

send  neutral  state  in  case  last  init 

clear  flag  from  IFR 

wait  for  neutral  state  to  xmit 

return  from  subroutine 
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/********************************************************************** 

/* 

/*  File:  HOSTAPPl.CPP  Source  code  for  host  application 
/* 

/********************************************************************* / 
#include  <HI54X.H> 

#include  <stdio.h> 

#include  <conio.h> 

#include  <stdlib.h> 

extern  int  datareg  []  ,  statreg  []  ,  ctrlreg  []  ; 
extern  int  pport,  portmode, Re adde lay; 


void  main (void) 

{ 

FILE  *fp; 

if  { {fp=fopen("datal.dat", "w") )==NDLL)  /*  Open  file  */ 

{ 

clrscr  0 ; 

printf ("Cannot  open  file  An"); 
exit (0)  ; 

} 

portmode=0;  /*  4 -bit  mode  */ 

Readdelay  =  20;  /*  m  case  host  slow*/ 

clrscr  0;  /*  dear  the  screen  */ 

if ( (pport =locate_port 0 )  >=  5) {  /*  Find  the  port.  */ 

printf("No  connection\n" ) ;  /*  If  no  connection  */ 

backoutO;  /*  then  leave  board  */ 

exit{0);}  /*  in  known  state  */ 

else{ } 

_setcursortype (_NOCDRSOR) ;  /*  Hide  text  cursor  */ 

set_latch(l,l) ;  /*  Keep  DSP  running  */ 

int  word  =0,  col=0;  /*  and  bring  PAL  out*/ 

/*  out  of  Tri-state  */ 

col=0; 

gotoxyd,!);  /*  go  to  home  */ 

send_word(0x0808,  C_SEND)  /*  dear  the  HINT  */ 

HINT(IOOOO);  /*  Wait  for  nxt  HINT*/ 

send_word{0xl200,  A_SEND) ;  /*  Goto  0x46  entries*/ 

/*  before  buffer  */ 
for (int  buf=0  ;  buf  <  0x271;  buf++) 

/*  change  here  if  fs  change*/ 

word  =  read_word (D_READ) ;  /*  Read  word  from  pp*/ 

printf ("%4.4x  ",  word);  /*  Print  it  to  scr  */ 
fprintf(fp,  "%d\n",  word);  /*  Output  to  file  */ 
if (col  >=  13) {  /*  in  14  columns  */ 

col=0; 

printf ("\n") ; } 
else{col++; } 

} 

_setcursortype  (_NORMALaiRSOR)  ;  /*  Ret  normal  cursor*/ 
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fclose (fp) ; 
backout ( ) ; 

exit (0) ; 


/*  Close  file  */ 


/*  Leave  board  in  */ 
/*  known  state  */ 


79 


/**********************************************************************/ 

/* 

/*  File:  H0STAPP2.CPP  Source  code  for  host  application 
/* 

/**********************************************************************^ 

# include  <HI54X.H> 

# include  <stdio.h> 
tinclude  <conio.h> 

#include  <stdlib.h> 

extern  int  datareg  []  ,  statreg  []  ,  ctrlreg  []  ; 
extern  int  pport,  portmode,Readdelay; 


void  main (void) 

{ 

FILE  *fp; 

if  { (fp=fopen{"data2.dat", "w") )==NDLL) 


/*  Open  file  */ 


clrscr 0 ; 

printf ("Cannot  open  file  .\n") ; 
exit (0) ; 


portmode=0;  /*  4-bit  mode  */ 

Readdelay  =20;  /*  In  case  host  slow*/ 

clrscr  0;  /*  clear  the  screen  */ 

if ( (pport=locate_port 0 )  >=  5) {/*  Find  the  port.  */ 

printf ("No  connect ion\n" ) ;  /*  If  no  connection  */ 


backout ( ) ; 
exit (0) ; } 

else{ } 

_setcursortype (_NOCURSOR) ; 
set_latch(l, 1) ; 
int  word  =0,  col=0; 


col=0; 

gotoxyd,!);  /*  go  to  home  */ 

send_word (0x0808,  C_SEND) ;  /*  Clear  the  HINT  */ 

HINT(IOOOO) ;  /*  Wait  for  nxt  HINT*/ 

send_word ( 0x1200 ,  A_SEND) ;  /*  Goto  0x46  entries*/ 

/*  before  buffer  */ 
for (int  buf=0  ;  buf  <  0x200;  buf++) 

/*  Change  here  if  fs  is  changed*/ 

word  =  read_word(D_READ) ;  /*  Read  word  from  pp*/ 

printf ("%4.4x  ",  word);  /*  Print  it  to  scr  */ 
fprintf(fp,  "%d\n",  word);  /*  Output  to  file  */ 
if (col  >=  13) {  /*  in  14  columns  */ 

col=0; 

printf ("\n") ; } 
else{col++; } 

} 


/*  then  leave  board  */ 
/*  in  known  state  */ 

/*  Hide  text  cursor  */ 
/*  Keep  DSP  running  */ 
/*  and  bring  PAL  out*/ 
/*  out  of  Tri- state  */ 
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} 


_setcursortype  {_NORMALCURSOR)  ; 
fclose (fp) ; 
backout ( ) ; 
exit  (0) ; 


/*  Ret  normal  cursor*/ 
/*  Close  file  */ 

/*  Leave  board  in  */ 
/*  known  state  */ 
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/***************************************************** 

/* 

/*  File:  MAIN.C  Source  code  for  main  program 
/* 

/***********************************************^**^^^^^^^^^^^^^^^^^^^^ 
#include  "math.h" 

#include  "stddef.h" 

#include  "stdlib.h" 

#include  "stdio.h" 

# include  "conio.h" 


/ 


tdefine  fsl  625 
#define  fs2  324 


main ( ) 

{ 

FILE  *fpl,  *fp2; 

int  dl,  d2,  xl[fsl],  x2 [f s2] ; 

int  k,  out,  freq; 

int  glde(int  k) ; 

double  xrol [fsl] ,  xiol [fsl] ,  xol [fsl] ; 

float  pi  =  3.1415926,  tpi; 

int  n,  u; 

int  i,  dftl,  nl; 

double  maxi ,  max2 ; 

double  xro2 [fs2] ,  xio2 [fs2] ,  xo2 [f s2] ; 
int  j,  dft2,  n2; 

if  ( (fpl=fopen( "data5.dat ", "rt") )==NDLL)  /*  Open  file  */ 

clrscrO  ; 

printf ("Cannot  open  file  .\n"); 
exit (0) ; 

} 

if  ( (fp2=fopen("data6 .dat" , "rt") ) ==NDLL)  /*  Open  file  */ 

clrscr  0 ; 

printf ("Cannot  open  file  .\n"); 
exit (0) ; 

} 

for (n=0;n<fsl;n++) 

{ 

fscanf(fpl,  "%e  ",  &dl)  ; 
xl [n] =dl ; 

} 

for (n=0 ;n<f s2 ;n++) 

{ 

fscanf(fp2,  "%e  ",  &d2)  ,- 
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x2 [n] =d2 ; 

} 

f close (fpl) ; 
f close (fp2) ; 

tpi=2*pi; 

for (u=0;u<fsl;u++) 

{ 

xrol [u] =0.0; 
xiol [u] =0.0; 
for (n=0;n<fsl;n++) 

{ 

/*--  Xr [u]  =  (1/fsl)  sum  {xr [n] .cos {2PI.u.n/fsl) }  --*/ 
xrol [u]  =  xrol [u]  +  xl [n] *cos (tpi*u*n/fsl) ; 
/*--  Xi [u]  =  -  (1/fsl)  sum  xr[n] .sin(2PI.u.n/fsl)  --*/ 
xiol [u]  =  xiol [u]  -  xl [n] *sin{tpi*u*n/fsl) ; 

} 

xrol [u] =xrol [u] /fsl; 
xiol [u] =xiol [u] /fsl; 

xol [u] =sgrt (xrol [u] *xrol [u] +xiol [u] *xiol [u] ) ; 

} 

dftl=0; 
nl=fsl/2+l; 
maxl=xol  [0] ; 
for  ( i=l ; i<nl ; i++ ) 

{ 

if (xol [i]  >  maxi) 

{ 

dftl=i; 
maxl=xol [i]  ; 

} 

} 

for (u=0;u<fs2;u++) 

{ 

xro2 [u] =0.0; 
xio2 [u] =0.0; 
for (n=0;n<fs2;n++) 

{ 

xro2 [u]  =  xro2 [u]  +  x2 [n] *cos (tpi*u*n/fs2) ; 
xio2  [u]  =  xio2 [u]  -  x2 [n] *sin(tpi*u*n/fs2) ; 

} 

xro2 [u] =xro2 [u] /fs2; 
xio2 [u] =xio2 [u] /fs2; 

xo2 [u]  =sqrt (xro2 [u] *xro2 [u] +xio2 [u] *xio2 [u] ) ; 

} 

dft2=0; 

n2=fs2/2; 
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max2=xo2  [0]; 
for ( j=l; j<n2; j++) 

{ 

if(xo2[j]  >  max2) 

{ 

dft2=j ; 
max2=xo2 [j] ; 

} 

} 

k=dft2-dftl; 

out=glde (k) ; 

freq=abs (dftl+out*fsl) ; 

printf("The  frequency  is  %d",  freq) ; 

for ( ; ; ) ; 

} 

glde (k) 

{ 

float  mult , bl ,  b2 ,  mtest ,  mdl ,  md2 ,  tnx; 

/*  This  section  solves  the  linear  diophantine  equation  fsl*bl  +  fs2*b2 
na  where  fsl  and  fs2  are  the  sampling  frequencies 

and  na  is  the  greatest  common  divisor  and  returns  the  value  bl,  b2  and 
na.  fsl  and  fs2  are  assumed  positive  */ 

float  bol,bo2,ma,na, irem,bo3,bo4; 
int  iquot ; 

bol=l; 
bo2=0 ; 
bl=0; 
b2=l; 

/*  Pl3ce  fsl  and  fs2  in  ma (dividend)  and  na  (divisor)  respectively  */ 

ma=fsl; 

na=fs2; 

/*  Calculate  quotient  and  remainder  */ 

iquot =ma/na; 
ir em=ma -na* iquot ; 

/*  If  remainder  is  not  zero,  reset  dividend  and  divisor  */ 

while  (irem>0) 

{ 

/*  calculate  new  coefficients  */ 

/*  redefine  bol,  bo2,  bl  and  b2  */ 


bo3=bol- iquot *bl ; 
bo4  =bo2  -  iquot  *b2  ,- 
bol=bl; 
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} 


bo2=b2 ; 

bl=bo3; 

b2=bo4 ; 

ma=na ; 

na=iretn; 

iquot=Tna/na; 

ireTn=ma-na*  iquot ; 


/*  redefine  dividend  and  divisor  */ 
/*  reapply  Euclidean  algorithm  */ 


/*  To  check  whether  the  equation  is  solvable,  na  must  be  a  factor  of  k 
for  the  equation  to  be  solvable.  */ 


mult=k/na; 
if  ( (k-mult*na) ==0) 
{ 

bl=bl*mult ; 

b2=b2*mult ; 
mtest=bl; 

mdl=fsl/na; 
md2=fs2/na; 
Tnx=bl  ; 
tnx=mx+md2 ; 


/*  Equation  is  solvable  */ 


/*  These  new  values  solve  the  */ 

/*  diophantine  equation  */ 

/*  To  check  whether  bl  and  b2  */ 

/*  are  the  least  values  that  */ 


/*  satisfies  the  diophantine  equation  */ 


return 

} 


while  ( (abs (mx) -abs (bl) )  <  0) 

{ 

bl=mx; 
b2=b2“mdl ; 

Tnx=tTix+md2 ; 

} 

if  ( (mtest-bl) ==0) 

{ 

mx=bl ; 
mx=mx“md2 ; 

while  ( (abs (mx) -abs (bl) ) <0) 

{ 

bl=mx; 
b2=b2+mdl; 
tnx=mx-md2 ; 


( {int)bl) ; 
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